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PHYS393 - Statistical Physics 
Part 3: The Maxwell- Boltzmann Gas 



The Boltzmann distribution 



In the previous parts of this course, we derived the Boltzmann 
distribution: 

N fl 

nj = —e kT, (1) 



We applied this distribution to two example systems: a spin-| 
magnetic solid, and a collection of 1-dimensional quantum 
harmonic oscillators. We saw how to derive a number of 
properties of these systems, including the energy, heat capacity, 
and entropy as functions of temperature. 



The partition function Z plays an important role in calculating 
quantities such as the total energy and entropy. It is defined by: 

Z = J>-3\ (2) 

j 
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The Maxwell-Boltzmann gas 



In this part of the lecture course, we shall apply the ideas and 
methods we have learned to a further system: the 
Maxwell-Boltzmann gas. This provides a simple model of a 
gas, that describes some important properties of gases quite 
well in a number of cases. 

Discussing this example will also reinforce some the ideas and 
methods we have learned so far. 

The most important new concept is that of the density of 
states. This is introduced because the energy levels in this 
system have a more complicated structure than in the previous 
two examples. We need to make some approximations to be 
able to evaluate the partition function and other important 
quantities. 
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The Maxwell-Boltzmann gas 



Consider a gas in a box of fixed volume. The gas consists of 
particles that are free to move within the box, but cannot 
penetrate the walls. We will assume that the density is low, so 
that collisions between gas particles are infrequent: this 
satisfies the condition that particles in the system interact only 
weakly with each other. 

Our goal, as in the previous examples, is to find quantities such 
as energy, heat capacity and entropy as functions of 
temperature. In principle, we use the same equations and 
techniques for this system as for the previous two examples. 
However, we first need to do some work to identify the energy 
levels. 

We shall use a quantum model for the gas particles: this will 
ensure that the allowed energy levels within the box occur at 
discrete intervals. 
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The Maxwell-Boltzmann gas 



Our system consists of a low-density gas of particles, which 
move freely apart from occasional collisions with the (rigid) 
walls or with each other. To distinguish this gas from others 
that we will consider later, we refer to the system as a 
Maxwell-Boltzmann gas. 

A further assumption we will make, at least for the time being, 
is that the gas particles have no internal structure. This means 
that there is no energy associated with rotations or internal 
vibrations of the particles. In practice, such a situation is 
closely approximated if the gas particles consist of single atoms 
- as in the case of noble gases (helium, neon...) 
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James Clerk Maxwell, 1831-1879 
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Maxwell-Boltzmann gas: wave functions 


For simplicity, consider a square box of side length L. The 




wave function ip(x,y,z) of a particle in the box must satisfy the 


boundary conditions: 




ijj(x,y, z) = 0 for x = 0, or y = 0, or z = 0, 


(3) 


and: 




ip(x,y,z) = 0 for x = L, or y = L, or z = L. 


(4) 


The boundary conditions are satisfied by a wave function of the 


form: 




?/>(x, y, z) = As\r\(k x x) s\r\(k y y) s\r\(k z z), 


(5) 


if: 

; 7rn x 1 irn y ; 7rn z 
— L ' — L ' — L ' 


(6) 


where n x , n y , and n z are positive integers. 
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Maxwell-Boltzmann gas: energy levels 


The wave function (5) is an energy eigenfunction, since it 




solves the time-independent Schrodinger equation: 






(7) 


where H is the Hamiltonian operator: 




H = V 2 

2m 


(8) 


and e is a constant, the energy of the particle. Substituting 


from (5), we find: 




m = L v2 * = L 


(9) 


Hence, the energy levels are given by: 




£ = 2mL 2 + n y + = kd + n y + ■ 


(10) 


The parameter 9 (with units of temperature) is defined: 




8- hV 
2mkL 2 


(11) 
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Maxwell- Boltzmann gas: energy levels 



We now run into trouble. If we want to evaluate populations, 
total energy, heat capacity etc. as functions of temperature, 
the next step is to evaluate the partition function: 

Z = J>-3\ (12) 

3 

Here, the index j will actually be a triple index, for %, n y and 
n z . But quite apart from that, the summation is now difficult 
to evaluate. 



To make things more convenient, we shall assume that we can 
replace the summation by an integral, and write: 

£ j roo £ 

Z = ^2e~kT g(£)e~kTd£. (13) 

Here, g(s) is the density of states: the number of energy states 
in the (small) energy range e to e + de is g{e)de. 
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Maxwell-Boltzmann gas: density of states 

To evaluate the partition function, we need an expression for 
the density of states, g(e). As an intermediate step, we will find 
the density of states in fc-space, g(k), and then convert to the 
density of states as a function of energy. 

In /e-space, states lie on a 3-dimensional rectangular grid of 
points with spacing ir/L between points. 

The number of points in the region bounded by k and k-\-dk is: 

, N A-7rk 2 dk 1 , N 

g(k)dkn— x— — 3. (14) 

8 (ir/L)° 

The first factor is the volume between spherical shells of radius 
k and k + dk, and the second factor is the number of grid 
points per unit volume of /c-space. The approximation is valid 
for large k, i.e. at high energies. 
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Maxwell-Boltzmann gas: density of states 




, 4-7rk 2 dk 
g(k) dk « x 



(tt/L)' 



(for large k). 
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Maxwell-Boltzmann gas: density of states 

With the assumption of large k, we have: 

4tt/c 2 



g(k) dk — 



and: 



Hence: 



8(tt/L) 3 
h 2 k 2 



dk, 



g(k) dk — 



2m 

A-7T 2me dk 



8(tt/L) 3 h 2 de 
Also, from equation (16), we have: 



de. 



k = 



I2me 



dk 
de 



1 2m 



2\h 2 t 



So, writing V = L 3 , we have for the density of states: 



4-7T 2me 1 2m V (2m\i 

g\e) de = — — „ 0 — \ 7o-de = — ~ —75- v e « e - 



8(tt/L) 3 ft 2 2V/i 2 e 4tt 2 V ft 2 



(15) 



(16) 



(17) 



(18) 



(19) 
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Maxwell-Boltzmann gas: partition function 

In terms of the parameter 0 (11), the density of states can be 
written: 

^=^5*' (20) 



we can evaluate the partition function (13): 

Z= / g(e)e~kTd£. (21) 
JO 

Substituting for g{e) de from (20) into (21) gives: 



7T 1 f°° £ 

Z = , / y/ee'Wrde. (22) 



The result of the integral is: 

J Vie'vrde = ^-(fcT)S. (23) 
So the partition function becomes, for high temperature: 

H# ^ 
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Maxwell-Boltzmann gas: low density approximation 

So far we have assumed that the gas is of low enough density 
that particles interact only occasionally. In a classical model, 
this simply means that the total volume occupied by the 
particles is much less than the total volume of the container. 

However, we are using a quantum model. In quantum 
mechanics, the strength of the interaction between two 
particles depends on how much their wave functions overlap. 
Two particles in states with similar wave functions will interact 
more strongly than the same particles in states with very 
different wave functions. 

Therefore, to validate our assumption of weak interactions, we 
need to check that the mean occupancy of any state in the 
system is low. 
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Maxwell-Boltzmann gas: low density approximation 



In the Boltzmann distribution, the number of particles in a 
state with a particular energy ej is: 

nj = —e kT. (25) 

Since the exponential factor can have a value up to 1, the 
occupancy of any state will be small if: 

-<1. (26) 

This is the condition that must be satisfied for our assumption 
of low density, and therefore infrequent collisions between 
particles, is to be satisfied. 
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Maxwell-Boltzmann gas: low density approximation 
Using equation (24) for the partition function, we can write: 

»= N (™y = »( 2 ^f. (27) 

Z \ttTJ V\mkT) 



Let us consider the case of (monatomic) helium gas, at s.t.p. 
One mole of an ideal gas at s.t.p. occupies 22.4 litres. 
Therefore, we use the following values: 

N = 6.02 x 10 23 , 

V = 22.4 x 10~ 3 m 3 , 

h = 1.05 x 10 -34 Js, 

m = 6.69 x 10~ 27 g, 

k = 1.38 x 10 -23 J/K, 

T = 273 K. 

With these values, we find N/Z « 3.9 x 10 -6 : the low density 
condition is satisfied. 
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Maxwell-Boltzmann gas: low density approximation 

Let us consider the worst case: helium gas at 1 atm pressure 
and 5 K (the boiling point is 4.2 K). In this case, the volume 
(assuming ideal gas behaviour) is: 

rpl c 

V< « —V = x 22.4 dm 3 . (28) 

T 273 

With the new values: 

V = 0.410 x 10~ 3 m 3 , 
T = 5 K. 

we find: 

— « 0.085. (29) 
Z V 1 

This is now getting somewhat close to 1: the 
Maxwell-Boltzmann theory should still give a reasonable 
description of the behaviour of helium at 5 K, but we should 
also not be surprised if the predictions of the model are not 
terribly accurate. 
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Maxwell-Boltzmann gas: total energy (1) 



Now that we have the partition function Z (for high 
temperature), we can calculate the total energy from: 

U = N—±-^, (30) 
dp 

where (3 = —1/kT. Writing the partition function in the form: 

Z=(-^)\ (31) 



we find: 



4k{30 



ln(Z) = 3 In (--— ). (32) 

v J 2 V 4kf36) v J 

Differentiating with respect to (3 gives: 

V = N d -^1 = -IE. (33) 

The result, valid in the high temperature (classical) limit, is: 

U = ^NkT. (34) 
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Maxwell-Boltzmann gas: total energy (2) 



Of course, we should also be able to calculate the total energy 
by summing the energies of the various states, weighted by the 
population determined from the Boltzmann distribution. 

The number of states within the energy range e to e + de is 
g{e)de. The number of particles within each of these states is: 

n(e) = — e «\ (35) 

and these particles have energy e. Hence, the total energy of 
the system is: 

N r°° s 
U = — g{e)ee kr de (36) 

2 / 7T \§ f°° 3 __e 



= — / eie'kTde. (37) 

7T 2 KkTj JO V } 

In the last line, we have substituted for Z and g(s) from 
equations (24) and (20) respectively. 
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Maxwell-Boltzmann gas: total energy (2) 



Now we need to perform the integral in equation (37). We 
make use of the standard result: 

/ x n e~adx = a n+1 r(n+ 1). (38) 
JO 

The gamma function l~(n + l) is a generalisation of the factorial 
function: for integer n, l~(n + 1) = n\. Using (38), we find: 

£ 2 e -feT d e = (kTp rl-j = (kT)2. (39) 

Substituting this result into (37) gives immediately: 

U = ^NkT, (40) 
in agreement with equation (34). 

Although we derived this expression for the energy assuming a 
square box, the result is valid for a box of any shape. 
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Maxwell-Boltzmann gas: total energy 



We made a high temperature approximation in order to 
evaluate the partition function: this limits our results to the 
classical regime. 



However, we can evaluate numerically the exact (quantum) 
expressions for the partition function and for the total energy. 
The total energy can be written in the form: 

,20 

(41) 



C/ = -E^e-^ = 3iV^ E ^" 2e ""' T 



z <-r 

3 



oo — n* 



Eoo 
n=l 



T 



Let us plot the exact results, and compare with the classical 
limit... 
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Maxwell-Boltzmann gas: total energy 




200 



400 600 

r/e 



800 1000 



At low temperature T < 6, we see the quantum effects as a 
non-zero energy of the particles in the ground state. 

At high temperature, T > 0, the energy approaches the 
classical limit U = ^NkT, shown by the dashed line. 
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Maxwell-Boltzmann gas: total energy 



At what temperature do we expect quantum effects to become 
important? Let us consider 1 mm 3 of helium gas. Using: 

m = 6.69 x lO -27 kg, 

L = 0.001m, 

h = 1.05 x 10 -34 Js, 

k = 1.38 x 10 -23 J/K, 



we find: 



fc2 2 

9 = ^ = 5.9 x 10" 13 K. (42) 

2mkL 2 V J 



This is an extremely low temperature, by anybody's standards. 
For all practical purposes, the classical theory should describe 
the behaviour of the gas very well. 
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Maxwell-Boltzmann gas: heat capacity 



From the expression (34) for the total energy of the gas: 

U = —NkT, 
2 

it is trivial to find the heat capacity (for T > 9): 

Cv=(P\ = 3 Nk. (43) 
v \dTJv 2 

The heat capacity at constant volume is constant. Note that in 
this case, it is important to apply the condition of constant 
volume, because the energy levels depend on the volume of the 
box. 



We should remember that these results are valid only in the 
high temperature limit. In this limit, the energy levels are close 
together compared to the typical energy of a particle (of order 
kT), so the gas behaves like a classical system. 
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Maxwell-Boltzmann gas: entropy 



Using the classical approximation for the partition function Z 
(24): 

Z=(^Y (44) 



40 

we can find an expression for the Helmholtz free energy: 

F = —NkT ln(Z), (45) 

and then find an expression for the entropy of the 
Maxwell- Botlzmann gas, using: 

'dF\ 



S = - 
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Maxwell-Boltzmann gas: entropy 



Using (24) we can write for the Helmholtz free energy: 

F = -NkT\n(Z) = -^NkT\n(^j - ^NkT\n(^j . (47) 



Differentiating with respect to T, we find that the entropy is 
given by: 

s =-(S.=l-[ i+in G)] + l™(D- ^ 



Again, this is an expression valid in the classical limit, T > 0; 
and, as for the case of the collection of harmonic oscillators, it 
has the form: 

S = S 0 + C v lnQn ( 49 ) 

that we expect from classical thermodynamics, for a system 
with constant heat capacity at constant volume. 
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Maxwell-Boltzmann gas: entropy 



Also note that, as in the case of the harmonic oscillators, the 
expression for the entropy includes Planck's constant, even in 
the classical limit. Entropy is inherently a quantum 
phenomenon. 

Comment: for an alternative discussion, using Boltzmann's 
equation S = k\n£2, see Guenault (second edition) page 69. 
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Maxwell-Boltzmann gas: further properties 



We can explore further properties of the Maxwell-Boltzmann 
gas. To complete our analysis of the monatomic gas, we shall 
calculate: 



• the pressure of the gas, and its relationship to temperature 
and volume; 



• the energy distribution of particles within the gas; 



• the velocity distribution of particles within the gas. 
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Maxwell-Boltzmann gas: pressure 



Pressure is a thermodynamic function, used in the description 
of macroscopic systems. We can derive the pressure from the 
Helmholtz free energy, in a similar way to the derivation of the 
entropy. 

Since the Helmholtz free energy F is given by: 

F = U -TS, (50) 
we can use the first law of thermondynamics in the form: 

dU = TdS - P dV, (51) 

to find that: 

dF = -SdT -pdV. (52) 

Hence, by considering a change in volume at constant 
temperature, we find: 

'dF\ 



V 
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Maxwell-Boltzmann gas: pressure 



Now we use: 

F = —NkT ln(Z). (54) 

For a monatomic Maxwell-Boltzmann gas, the partition 
function Z is given (in the classical limit) by (24): 



Using (11): 



we can write: 



Z=(^) ■ (55) 



0 = (56) 



2mkL 2 ' 



mkTL 2 \ 2 V 



z = lwJ =v (57) 



where: 
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Maxwell-Boltzmann gas: pressure 



Since we have: 



V 

Z = ~, (59) 



we find: 

P 



\dv) T \ dV ) T v K J 



Thus, the pressure, volume and temperature for a monatomic 
Maxwell-Boltzmann gas are related by: 

P V = NkT. (61) 
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Maxwell-Boltzmann gas: pressure 

The equation of state (61) for a monatomic 
Maxwell-Boltzmann gas: 

P V = NkT, 

corresponds to the ideal gas equation: 

pV = nRT. (62) 

Hence, a monatomic Maxwell-Boltzmann gas behaves (in the 
classical limit of high temperature) as an ideal gas. 

Furthermore, since N is the number of particles in the gas, and 
n the number of moles, N/n = A, where A is Avogadro's 
constant. Thus, we have a relationship between Boltzmann's 
constant k, and the gas constant R: 

R = kA. (63) 
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Maxwell-Boltzmann gas: energy distribution 



Previously, we found the density of states in a 
Maxwell-Boltzmann gas (20) as a function of energy: 



g(e) de = 



7T 



4 V (k6) 3 



de, 



(64) 



We also have the partition function, Z (24): 



\46) 



(65) 



Hence, using the Boltzmann distribution, the number of 
particles n(e)de in the energy range e to e + de is: 



N ^_ e 

<e) de = -g(e)e W de = Jj-^ 



e kr de. 



(66) 
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Maxwell-Boltzmann gas: energy distribution 




s/kT 
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Maxwell-Boltzmann gas: energy distribution 
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Maxwell-Boltzmann gas: energy distribution 



The peak in the distribution occurs at: 

dn(e) 



de 



= 0, 



which we find gives: 



e = -kT. 
2 



(67) 



(68) 



The average energy over all particles is given by: 

1 f°° 



— — / en(e) de. 

nJo 



(69) 



We have, of course, evaluated this integral before, when 
calculating the total energy. The result, using equation (38), is: 



e = -kT. 
2 



(70) 
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Maxwell-Boltzmann gas: velocity distribution 

Finally, we turn to the velocity distribution of particles in a 
monatomic Maxwell-Boltzmann gas. We can obtain this 
directly from the energy distribution, by substituting: 

1 o 



e = -mi) , 
2 



de = mv dv 



and: 



into equation (66). 



The resulting expression for the energy distribution is: 



(71) 



(72) 



n(v) dv = A-7rN 



m 



2<nkT 



2 q mv 

V^e 2kT dv. 



(73) 



The velocity distribution n(v) is known as the 
Maxwell-Boltzmann distribution. 
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Maxwell-Boltzmann gas: velocity distribution 



T 1 <T 2 < T 3 
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Maxwell-Boltzmann gas: most probable velocity 



The most probable velocity for a particle in the gas is the 
velocity at which the velocity distribution peaks. This is found 
by solving: 

dn 

— = 0. (74) 

dv 



Since: 



we find: 



3 2 

mv 



n(v) = 4tt (^^) 2 ^ 2 e 2fcr , (75) 



— = 0 when v = v = \ . (76) 

dv V m 
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Maxwell-Boltzmann gas: mean velocity 



The mean velocity of particles in a Maxwell-Boltzmann gas is 
given by: 

1 fOC f m \ 2 f°° -a _rmPi 

v = — v n(v) dv = 47r / re 2kT dv. (77) 

NJo W \27rkTJ Jo V y 



The integral is the third Maxwell-Boltzmann integral, I3 (see 
Appendix A). It is evaluated using: 

roo „ , 2 1 1 

I 3 = / x 3 e- 6x = /i = -A*. (78) 



Hence, putting b = m/2kT, we find: 



-— ( m Y 1 ( 2kT \ 2 ~ l 8kT 

V ~ 77 \27rkTJ 2V~m~/ V™~' 
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Maxwell-Boltzmann gas: mean square velocity 

The mean square velocity of particles in a Maxwell-Boltzmann 
gas is given by: 

v 2 = v z n(v) dv = 4tt / e 2fcT dv. (80) 

7V7o W \27rkTJ Jo V y 



The integral is the fourth Maxwell-Boltzmann integral, I4 (see 
Appendix A). It is evaluated using: 

A . 3 



/ 4 = x 4 e~ bx2 dx = ^I 2 , (81) 



and: 



2 7o 26 U 46V 5 V y 



Hence, 



7 4 = -"L/^, (83) 
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Maxwell-Boltzmann gas: mean square velocity 



Thus, putting b = m/2kT, we find: 



^ 4l m^al^. (84) 

\2irkTj 8 V m J m K J 



Note that this implies that the mean kinetic energy of a 
particle in the gas is: 

1 3 

i = -mv 2 = -fcT, (85) 

which is consistent with our result for the total energy of the 
gas (34): 

U = ^NkT. (86) 
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Maxwell-Boltzmann gas: velocity distribution 



In terms of the parameter: 



IkT 



m 



(87) 



the most probable velocity v, mean velocity v and root mean 
square velocity \ftP- are given by: 



v = V2vt, 



(88) 



'8 

v = \l —vj>, 

IT 



— v3vt- 



(89) 
(90) 
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Maxwell-Boltzmann gas: velocity distribution 
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Diatomic Maxwell-Boltzmann gas 



So far, we have assumed that the particles in the gas have no 
internal structure. In other words, the only kind of energy they 
can possess is the "translational" kinetic energy associated 
with their velocity as they move through the box. 

This is a real model for monatomic gases: but there are only a 
few examples of such gases. Most gases are composed of 
molecules consisting of two or more atoms. 

The internal structure in a diatomic molecule provides 
additional energy states for the molecule. These energy states 
are associated with rotation and vibration of the molecule. We 
also now include the possibility of electronic excitation of one 
or both of the atoms in the molecule. 
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Diatomic molecule degrees of freedom 




rotation electronic excitation 
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Diatomic Maxwell-Boltzmann gas 



Note that rotational energy states only include those with 
orbital angular momentum, in which the positions of the atoms 
change with time. Generally, the intrinsic angular momentum 
associated with the spin of the atoms is fixed. 




Orbital angular momentum. 
Energy levels: s, = ^y / ( / + 1 )- 




Spin angular momentum. 
Fixed energy. 
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Diatomic Maxwell-Boltzmann gas 



We must modify our theory to take into account the extra 
energy levels associated with rotations, vibrations, and 
electronic excitation. To keep things simple, we will assume 
that the degrees of freedom are independent of each other. 
That is, the vibration of the atoms within a molecule does not 
affect the rotational energy, and vice-versa. 



In practice, this is not exactly true, but is a good 
approximation. It allows us to write the total energy of a 
particle in the gas as a simple sum of the energies associated 
with the various degrees of freedom: 



£ — £ trans + £ rotn + £ vibn + e elec- 



(91) 
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Diatomic Maxwell-Boltzmann gas: partition function 

The partition function Z then becomes: 

Z = £ exp = £ exp (-^ ( et + £r + £v + e e )) • (92) 

The summations extend over all possible states. 

The partition function can now be expressed as: 

z = J2 • £ • £ • £ e_ ^- ( 93 ) 

Hence, the partition function can be written as a product of 
partition functions associated with each of the degrees of 
freedom: 

Z = Ztrans ' ^rotn ' ^vibn ' ^e\ec- (94) 
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Diatomic Maxwell-Boltzmann gas: partition function 

Recall that the energy and the entropy depend on the 
logarithm of the partition function: 

d\nZ 1 , N 

U = N , (3 = ; (95) 

d/3 ' H kT v J 

and: 

'dF\ 



S=- 



( — 1 , F = -NkTlnZ. (96) 

\dTJv 



Since Z is the product of partition functions associated with 
each of the degrees of freedom, the logarithm of Z becomes a 
sum of the logarithms of the different partition functions: 

InZ = lnZ t + lnZ r + In Z v + lnZ e . (97) 



This is very convenient. It means that the total energy, the 
entropy, and the heat capacity are simply sums of contributions 
from the different degrees of freedom. 
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Diatomic Maxwell-Boltzmann gas: heat capacity 

The heat capacity is something that is easy to measure directly 
in the laboratory; so we shall focus our attention on this. Our 
aim is to calculate the total heat capacity of a diatomic gas, 
taking into account all the degrees of freedom of the particles. 

Let us consider the contribution to the heat capacity of each of 
the degrees of freedom in turn, starting with the translational 
energy. This is very easy: everything is the same as in the case 
of the monatomic gas. 

The contribution to the heat capacity from the translational 
energy is: 



C$p = -Nk. (98) 



<V 2 

Remember that this is valid in the classical regime, which is 
appropriate for all gases under ordinary conditions. 
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Next, let us consider the electronic excitation. Typically, the 
first excited state will have an energy that is of order 10 eV 
above the ground state. The population of the first excited 
state, relative to the ground state, will be: 

e kT gj-eO 

~Z^ = e kT . (99) 

e kT 

where e\ - eq « 10 eV. At room temperature, T « 290 K, and 
kT « 0.025 eV. Hence, the first excited state has a population 
of approximately e -400 « 0 relative to the ground state. 

At ordinary temperatures, the electrons in all the atoms will be 
in the ground state: we expect electronic excitation to make no 
contribution to the heat capacity of the gas. 
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Now let us consider the vibrational motion. For this degree of 
freedom, we can treat the diatomic molecule as a quantum 
harmonic oscillator, with energy levels: 

ej =jhuj, (100) 

where u is determined by the strength of the bond between the 
atoms, and the mass of the atoms. 

We define a parameter 0 V , with units of temperature, such that: 

k6 v = hw. (101) 

From Part 2 of this course, where we considered a collection of 
quantum harmonic oscillators, we can then write for the 
vibrational contribution to the heat capacity: 

4 W = 4 ,2 - (102) 
eT - 1 
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Studying the behaviour of the heat capacity is one way to 
determine the value of the parameter 0 V for particular gases: 
this will become clear shortly, when we sum the contributions 
to the heat capacity from the various degrees of freedom. For 
now, we simply note that for many gases, the value of 0 V is of 
order of a few thousand kelvin: 



Gas 


0v 




o 2 


2000 


K 


N 2 


3000 


K 


H 2 


6000 


K 



Recall that for T <c 0v, the first excited energy level of the 
oscillator is significantly above the energy kT\ hence, at room 
temperature where T« 290 K, the molecules of most diatomic 
gases will be in their ground vibrational state. At room 
temperature, the vibrational degree of freedom makes little 
contribution to the heat capacity. 
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For T > 0 V , the heat capacity of a collection of quantum 
oscillators approaches the value Nk. However, at a few 
thousand kelvin, the molecules of many gases start to 
dissociate into individual atoms. This means that the 
vibrational degree of freedom only makes a contribution to the 
heat capacity over a relatively small range of temperature. 
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Finally, let us consider the rotational degree of freedom. From 
quantum theory, we know that the energy levels of a system 
with orbital angular momentum are given by: 

h 2 

ei = —t(t+l), (103) 

where I is a positive integer {t = 1,2,3...), and / is the moment 
of inertia. 

We also know that within each energy level, there are 

g = (21 + 1) distinct states (in other words, the energy levels 

have a degeneracy 2t + 1.) Thus, the partition function can be 

written: 

Z x = J>e~^ = £(2^ + 1) e -^ £+1) T. (104) 

t I 

where we define the parameter 0 r (with units of temperature) 
such that: 

h 2 

k6 r = . (105) 
21 V y 
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As usual, the parameter Or characterises the temperature at 
which the degree of freedom starts to make a contribution to 
the heat capacity. It depends on the moment of inertia of the 
molecule: this may be estimated from a classical model, if the 
masses mi and 777,2 and distance 2r between the atoms are 
known: 

J= mim2 r 2 . (106) 

777-1 + m 2 

Or, 6 r may be obtained from experiment. We will not go into 
the details, but simply quote some results: 



Gas 


Or 




0 2 


2.1 


K 


N 2 


2.9 


K 


H 2 


85 


K 


HD 


64 


K 
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The values of 6> r are very low: typically (except for hydrogen) 
below the boiling point of the gas. This is a consequence of 
the fact that rotational energies are quite low, certainly 
compared to vibrational energies. Typically, the rotational 
degree of freedom will always contribute to the energy of a gas, 
and hence to the heat capacity. 

Since at room temperature, most gases have T > 0 r , a gas at 
room temperature will be in a classical regime for the 
rotational energy, where the difference between energy levels is 
small compared to the typical rotational energy of a molecule. 
In other words, we can treat the energy levels as a continuum. 
This is fortunate, since otherwise it would be very hard to do 
the summation in the partition function... 
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Making a classical approximation, we can write the partition 
function for the rotational degree of freedom (104): 

Z r = J^(2£ + 1) e -^+l)T « ^°°(2^ + 1) e-^+^T (107) 



The integral is easier than it looks. We define a new variable y. 

y = £(e+l) .-. dy=(2l+l)dl. (108) 

Then the integral becomes: 

roc 9 r X 
Z r « / e - y Tdy = -. (109) 
JO ur 
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Now that we have the partition function (in the classical limit), 
we can find the total energy U r associated with the rotational 
degree of freedom. This is given by: 

Ur = N dU ^, f>=-± (HQ) 
Using Z x = ^ we find: 

Ur = NkT, (111) 
and hence the heat capacity is given by: 

C^p = Nk. (112) 

It is important to remember that this expression is valid in the 
classical limit, T > 0 r . For T <c 0 r , the rotational degree of 
freedom will not be excited, and the total energy and the heat 
capacity will be zero. 
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In summary, we expect the heat capacity of a diatomic gas to 
vary with temperature as follows: 



T range 








Total CV 


6 < T < 9r 




0 


0 




Or < T < 6y 


§iVfc 




0 




6y<T 


§iVfc 


7V/c 


Nk 


^Nk 



Recall that: 



2_2 



6> = 



H Z 7T 



2 ' 



6v = 



fojj 



(113) 



2/fc' * k ' 

where V is the volume of the container; / the moment of 
inertia of a gas molecule; and u> the frequency of vibration of 
the atoms in the molecule. 
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For real gases, the theory appears to work reasonably well.. 



4.5 




6000 
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The heat capacity of hydrogen 



A closer look at the heat capacity of hydrogen suggests that 
there may be a problem. A careful calculation suggests that, 
for the rotational heat capacity C^p , we should find curve A in 
the figure below; in fact, experimentally, curve C is found. 




T/a 
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The heat capacity of hydrogen 



The reason for the discrepancy is to do with the spins of the 
nuclei of the atoms in a hydrogen molecule. The nucleus of a 
hydrogen atom is a proton, which has spin-^. For a hydrogen 
molecule, the spins can be parallel, in which case they sum to 
spin-1 (ortho-hydrogen); or the spins can be antiparallel, in 
which case they sum to spin-0 (para-hydrogen). 

The spin-1 states have an additional three-fold degeneracy 
compared to the spin-0 states. We need to take this into 
account when finding the partition function. 
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The heat capacity of hydrogen 



In a hydrogen molecule, the total spin and the orbital angular 
momentum of the protons couple in an interesting way. 
Because protons are fermions, the wave function of the system 
must be antisymmetric under interchange of the coordinates of 
the two protons, i.e. 

ip(x 1 ,x 2 ) = -</>(x2,xi). ( 114 ) 

A state with orbital angular momentum £ is symmetric under 
interchange of the protons if £ is even, and antisymmetric if £ is 
odd. The state with spin-0 (spins of the two protons 
antiparallel) is antisymmetric; while the three states with spin-1 
(spins parallel) are symmetric. 

Thus, the possible combinations of spin and orbital angular 
momentum are: 

• £ even; spin-0 (single spin state for each £)\ 

• £ odd; spin-1 (three spin states for each £). 
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Taking into account the additional degeneracy associated with 
the spins of the hydrogen nuclei, the partition function for the 
rotational degree of freedom should become: 

Zv = £ (2£+l)e-^+ 1 )? + 3 £ {2£ + l)e~^ +1 ^ . (115) 
I even I odd 

If we compute the rotational heat capacity using this new 

partition function... 
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...we get curve B in the plot. There is an even great 
discrepancy with the experimental results, curve C! 



o 



B \ 






i 



TM 
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The heat capacity of hydrogen 



The reason for the discrepancy is that in hydrogen gas, the 
probability for a molecule to change from ortho-hydrogen to 
para-hydrogen (or vice-versa) in a collision with another 
molecule is extremely low. Molecules tend to stay in their spin 
states. If we take that into account, then we should treat the 
two components as different gases, each with their own 
rotational heat capacity, derived from the two distinct 
summations in the partition function (115). 

If ortho-hydrogen has three-times greater concentration in a 
sample than para-hydrogen (as may be assumed from there 
being three times as many possible states for any given 
rotational energy), then the total heat capacity will be: 



r (r) _ 3 r (r) , ! r (r) 



(116) 



This heat capacity gives the curve C in the previous figure... in 
agreement with experiment, at last! 
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The heat capacity of hydrogen 



Finally, we mention that if the measurement of the heat 
capacity is done slowly, or in the presence of a catalyst that 
increases the rate of conversion between the ortho- and para- 
states, then curve B in the previous figure is indeed obtained 
experimentally. This is because the system is able to achieve an 
equilibrium between the states, in which case the assumptions 
leading to the partition function (115) are satisfied. 

It is worth noting that the strange behaviour of the heat 
capacity of hydrogen really is a consequence of the fact that 
the atomic nuclei are idential fermions. A gas consisting of 
molecules of HD (an atom of hydrogen bound to an atom of 
deuterium) shows none of the strange behaviour of H 2 ; the 
variation of heat capacity of HD with temperature is similar to 
other diatomic gases. 
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Summary 



You should be able to: 

• Explain what is meant by a Maxwell-Boltzmann gas. 

• Derive expressions for the density of states and the partition function in 
a Maxwell-Boltzmann gas, in the classical limit. 

• Explain and check in particular cases the low-density condition for a 
Maxwell-Boltzmann gas. 

• Derive, as functions of temperature, the energy; heat capacity; entropy; 
and pressure of a Maxwell-Boltzmann gas. 

• Derive the energy and velocity distributions, as functions of 
temperature, for a Maxwell-Boltzmann gas. 

• Calculate the most probable velocity, mean velocity, and mean square 
velocity of particles in a Maxwell-Boltzmann gas. 

• Describe the additional degrees of freedom in a diatomic gas compared 
to a monatomic gas, and explain how they contribute to the partition 
function, and to the heat capacity. 

• Give a qualitative description of the behaviour of the rotational heat 
capacity of hydrogen as a function of temperature, and explain this 
behaviour in terms of the ortho and para states of hydrogen molecules. 
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To find the mean velocity and mean square velocity of particles 
in a Maxwell-Boltzmann gas, we need to evaluate integrals of 
the form: 



i n = r x u e - ux ~dx. 

Jo 

The case n = 0 is the standard Gaussian integral: 



Iq = e ox dx - 



2 V 6 



(117) 



(118) 



2 2 

For n=l, the integrand xe~ bx is the derivative of — e~ bx /2b, 



and so: 



oo 



I\ = / x e hxl dx = 
'0 



2b 



oo 



0 



1 

26' 



(119) 



Other integrals in the series may be found from a recursion 
formula, which is obtained by integrating by parts... 
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Appendix A: Maxwell-Boltzmann Integrals 



Integrating by parts gives: 

„n+l t°° 



In = \^ x n e~ hxl dx = 



JO 



x ( .-6.c 2 



71+ 1 



0 .0 72+ 1 



(-2bx)e- bx2 dx. 



Hence we find: 



. 00 2b n+2 _ bx 2 
I n = / — — x n ^ e ox dx = 



26 



'0 n + 1 
This may be written as 



n+ 1 



I n-\-2 — 



72+1 

26 



or, substituting n - 2 for n: 



In -^b~ In ~ 2 - 



(120) 

fn+2- ( 121 ) 
(122) 

(123) 
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